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Abstract. We investigate the first and second moments of the inverse 
participation ratio (IPR) for all eigenvectors of the Laplacian on finite random 
regular graphs with n vertices and degree 2 :. By exactly diagonalizing a large 
set of 2 :-regular graphs, we find that as n becomes large, the mean of the 
inverse participation ratio on each graph, when averaged over a large ensemble of 
graphs, approaches the numerical value 3. This universal number is understood 
as the large-n limit of the average of the quartic polynomial corresponding to 
the IPR over an appropriate (n — 2)-dimensional hypersphere of The graph 
ensemble mean variance of the inverse participation ratio for all graph Laplacian 
eigenvectors deviates from its continuous hypersphere average due to large graph- 
to-graph fluctuations that arise from the existence of localized modes. 
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Much of condensed matter physics involves the study of either localized or itinerant 
degrees of freedom that exist on the sites of a Euclidean lattice, defined by a notion 
of physical distance between a site and some number of proximate or “neighboring” 
ones. The relationship between a site and its neighbors defines both the dimension 
of space, d and a finite set of lattice vectors aj G j = 1,... ,d that can be used 
in combination with a set of integers {m} with G 1/“ to index n lattice sites via 
R. = Examples include the fourteen Bravais lattices in three spatial 

dimensions. However, it is often instructive to consider the same physical degrees of 
freedom on a non-Euclidean lattice, or graph, where no distance metric exists. A finite 
lattice of n sites is replaced by a graph G, consisting of a set of vertices V = {vi}, 
each connected by undirected edges to its neighboring vertices. The quantity Zi is 
known as the degree of the vertex Vi, and is equivalent to twice the spatial dimension 
for a hypercubic lattice. 

Studying a given physical model on a graph offers many technical benefits 
including the ability to: (i) study arbitrarily long range interactions where exact mean- 
field solutions may be available, (ii) smoothly tune and control the local dimensionality 
and (iii) easily encode the randomness and disorder that often exists in real systems. 
Celebrated examples from statistical physics include the solution of the Ising model 
and its generalizations; graph coloring and random percolation problems (for a review 
see Ref. mi Anderson’s model of non-interacting electrons hopping on a disordered 
lattice [2] was first solved on the Cayley tree (Bethe lattice) [3], providing deep physical 
insights into the nature of localization in quantum mechanical systems. More recently, 
the ability to study graphs in the limit z —>■ oo has lead to the development of 
Dynamical Mean Field Theory [3], allowing for systematic investigations of candidate 
microscopic models of the high temperature superconductors [S] . 

The discrete Laplacian matrix L plays a crucial role in defining any physical model 
on a graph, as it quantifies the kinetic energy cost of rapidly varying some degree of 
freedom amongst a set of neighboring vertices. For example, it encodes the classical 
dynamics in random vibrational networks [B] as well the onset of ferromagnetism in 
the classical [7] and quantum 0(n) model 0 E] on graphs. It appears in models 
of non-interacting bosons hopping between graph vertices, where the existence of a 
Bose-Einstein condensation transition on complex networks can be rigorously proven 

m- 

In this paper, we are interested in the properties of the Laplacian defined on 
finite sized regular graphs, where each of the n vertices is connected to exactly z 
neighbors. Examples with n = 100 and z = 3,8 are shown in Fig. Such graphs 

possess many important mathematical properties [12] while still retaining similarities 
to the physically realizable Bravais lattices discussed above. 

Much is known about the spectrum of the Laplacian on regular graphs, both 
in the thermodynamic limit n —> oo min] as well as on finite sized graphs un¬ 
substantially less is known about the eigenvectors of L [n na ng nn , as it is 
not possible to apply many of the standard tools of analysis for Euclidean lattices, 
including the Fourier transform. To address this gap, we systematically study the 
eigenvectors of the Laplacian matrix on a large class of finite size random regular 
graphs through brute-force numerical diagonalization. By investigating the statistics 
of the inverse participation ratio, we numerically observe that its mean across all 
eigenvectors approaches a finite universal value equal to 3, independent of graph 
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z = 3 
n = 100 



z = 8 

, n = 100 


Figure 1. Regular graphs with n = 100 vertices and degrees 2 = 3 and z = 8 
constructed via the Steger-Wormald algorithm [Ill- 


degree. We quantify the second moment of the distribution and find a dependence 
both on degree z and the number of vertices n. 

The paper is organized as follows: we begin with a formal definition of the discrete 
Laplacian on graphs in Section and describe our numerical results for the inverse 
participation ratio in Section [3| In Section we analyze the inverse participation 
ratio as a polynomial function on an (n — 2)-dimensional subsphere of n-dimensional 
real space. This perspective allows us to calculate exact values for the first and 
second moments of the inverse participation ratio over a continuous domain which 
contains the (terminal points of the) eigenvectors of the Laplacian. Sectionj^compares 
the values of the theoretically derived moments to those numerically computed from 
a large set of finite size random regular graphs. We analyze the deviations from 
the theoretically derived second moment with a discussion of localized eigenvectors 
and highlight implications for their use in computing physical observables on random 
regular graphs. 

2. The Graph Laplacian 

The Laplacian matrix generalizes the continuous Laplace operator A = V • V to 
encode variations of any continuous function (/> : V —> C which can take a value (fi 
on the vertex Vi. The physical importance of this matrix stems from the fact that 
solutions of A(() = 0 correspond to the Dirichlet energy functional which is stationary 
in some spatial region. The particular extension of A to a graph that we employ 
arises from the discrete approximation to the second continuous derivative of </) on a 
hypercubic lattice in d spatial dimensions with unit lattice spacing: 


d 



( 1 ) 


where e^- are the Cartesian unit vectors with elements e^j = Skj and S^j is the 
Kronecker i5-function. On a regular graph G consisting of n vertices Vi, each with 
degree z, the local connectivity is encoded in an adjacency matrix Aij where 



1 ; if Ui and Vj share an edge, 

0 ; otherwise. 


( 2 ) 
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Comparing with Eq. Q, we can write the elements of the graph Laplacian matrix as 
the difference between the degree and adjacency matrices of G: 

Lij = z6ij ^ij ( 11 ) 

and observe that z corresponds to twice the dimension on a hypercubic lattice. As 
mentioned in the introduction, L may appear in the Hamiltonian of numerous physical 
systems defined on a graph in the form: 


1 , 

H = (4) 

A spectral decomposition of L provides a route to the determination of the equations 
of motion governing a classical system, or the nature of the wavefunctions and allowed 
energy eigenstates for a quantum mechanical one. 


2.1. Exact diagonalization 

An ensemble of random regular graphs with with n vertices and degree z can be 
generated using the O(nz^) algorithm of Steger and Wormald [TT]. From the vertex 
neighbor list of each graph, we construct the adjacency matrix A and then exactly 
diagonalize the resulting n x n sparse Laplacian matrix L given in Eq. (§. In this 
paper, we present results for graphs with 

z e {3,4,5,10,15,20,25,30,35,40,45,50} 
n G {200,300,400,500,1000,2000,3000,4000,5000,10000} 


where all averages are performed over a set containing Ng graph realizations, with 
Nq = 5000 for n < 5000 and Nq = 1000 graphs for n > 5000. The exact number 
of unique graphs, Afc with a given n and z grows quickly with n but is unknown in 
general. An asymptotic result for degrees satisfying z < •y21ogn — 1 was proved by 
Bollobas in Ref. [5D]. Explicit counts for small n and z can be found at the Online 
Encyclopedia of Integer Sequences |2T], e.g. for z = 3 and n = 16, Mg = 4060. 

We begin our analysis by describing the eigenvalue distribution of L. For n >> 1, 
the limiting form of the density of states p(e), the probability of an eigenvalue falling 
between £ and £ + de, is given by the Kesten-McKay law mm]: 


P(e) 


-J(e) + 
n 


z •\/4(z — 1) — (e — z)2 
27r z^ — (e — z)2 


|e — z| < 2\/z — 1. 


( 5 ) 


For finite values of n, Metz et al. have recently computed the 1/n corrections to this 
expression, originating from the contributions of loops of all possible lengths. 

For a graph with n vertices, the eigenvalues of the Laplacian matrix {si} are 
determined by exact diagonalization, and a comparison to Eq. ([^ can be made by 
numerically constructing the histogram: 


(p(e)) 


71 ^ 


2=1 



( 6 ) 


where an average over Ng graph realizations is indicated by the angle brackets: 
{■■■) = (1/A^g) X)g(’■ ■ )■ results for n = 1000, z = 3,4,5,10,20,50 and 
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Figure 2. The probability of graph Laplacian eigenvalues for n = 1000 
vertices with degrees 2 : = 3, 4, 5,10, 20, 50 computed by diagonalizing numerically 
constructed graphs (shaded region) compared with the large-n limit Kesten- 
McKay law defined in Eq. (solid line). Eigenvalues in the continuum are 
bounded between Smin = z — 2\/z — 1 and Smax = 2 + z — 1, while the Perron- 
Frobenius mode with weight 1/n is shown as a spike at £ = 0. 


Nq = 5000 are shown in Fig. We observe only small graph-to-graph variations 
and find excellent agreement with the Kesten-McKay semi-circle law of Eq. ([^ using 
50 eigenvalue bins (solid line). For n > 1000 there is no visible discrepancy on this scale 
between the numerical results and the prediction for the large n limit. For finite sized 
random regular graphs, the spectrum of the Laplacian consists of a single eigenvalue 
at £ = 0 separated by a z-dependent gap [55] to a quasi-continuum of eigenvalues 
bounded between emin = z — 2^Jz — 1 and £max = z + 2yJz — 1. 

3. The Inverse Participation Ratio 

Less is understood about the set of eigenvectors E = {x : Lx = ex} HZlllHllllj 
although for certain classes of regular graphs with z ^ 0(n), they are believed 
(with high probability |23j 1 to be delocalized - meaning they have many non-zero 
components. The eigenvalue e = 0 with weight 1/n in Eq. ([^ and Fig. [^corresponds 
to the special case of the Perron-Frobenius mode: p = (l/-y/n,..., 1/yn) and via 
orthogonality it follows that x ■ p = y) ■ Xi = 0 for any eigenvector x ^ p. We wish to 
develop an understanding of the properties of the remaining eigenvectors £ = i?\ {p}, 
and in particular, determine how the non-zero elements of an arbitrary x G £ are 
distributed amongst its n coordinates. 

To this end, we study the notion of localization of an eigenvector using the inverse 
participation ratio. Historically, the participation ratio p was introduced to aid in 
classifying the properties of atomic vibrations in disordered lattices (53). It describes 
the fraction of the total number of sites which participate in a given normal mode 
vibration corresponding to the eigenvector x = {xi,..., Xn) € K" and takes the value 


P = 




( 7 ) 
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where /i’’ = X]r=i t)e thought of as the r**' moment of the kinetic energy 

of the mode. If a given normal mode only involves the motion of a single atom, it is 
characterized as localized and has p = 1/n. A vibrational mode consisting of all atoms 
participating equally is called extended and has p = 1. An equivalent measure was 
employed by Visscher [5S] to study the degree of localization of electronic eigenstates 
in the Anderson model [2] with implications for the existence of a metal-insulator 
(delocalization-localization) transition in the presence of disorder. 

When considering normalized eigenvectors x ■ x = ||a;|| = 1, it is often more 
convenient to consider the associated inverse participation ratio (IPR): 

n 

p-'^ = IPR{x) = (8) 

i^l 

For the Laplacian matrix in Eq. ([^ with x £ E we have 

1 < IPRix) < ^ (9) 

with bounds corresponding to the extended Perron-Frobenius (lower bound) and 
localized (upper bound) modes, respectively. The finite size scaling of IPR[x)/n 
as n —>■ oo provides information on the existence of a mobility edge, which defines 
the portion of the spectrum with robust delocalized states and has been extensively 
studied for a large class of random matrices [IHl H?! HH] . 

The inverse participation ratio thus provides a convenient single scalar value 
measuring the degree to which a particular eigenvector is localized (p“^ ~ 0{n)) 
or extended {p~^ 0(1))- To obtain information on the reduced set of Laplacian 

eigenvectors £ corresponding to non-zero eigenvalues, we construct a histogram of 
values in analogy with Eq. §. The non-linear form of the IPR necessitates that the 
order of averaging is important, and we must compute 

p(IPR) = 5(IPR - IPR(x)) (10) 

for each graph realization separately before averaging over graphs. The resulting 
graph averaged distributions are shown in Fig. for n = 1000, z = 3,4, 5,10, 20, 50 
and No = 5000. The solid lines show the results of a fit to a Gaussian distribution 
for each n and z and there are clear deviations which skew to larger IPR values, most 
notably for small z. For fixed n, increasing z decreases the width of the distribution, 
and slightly improves the residual of the Gaussian fit, but quite strikingly, the mean 
stays at a value that is numerically very close to 3, with the result appearing to be 
exact as n —>■ oo. This empirical finding warrants an explanation, which we provide 
by moving away from the discrete graph setting and employing tools of continuous 
analysis. 

4. Analysis of the IPR on a Hypersphere 

In this section, we take advantage of the geometry of the set of eigenvectors £ = 
E \ {p}. Their terminal points lie on a hypersphere S{p) which is orthogonal to 
the Perron-Frobenius mode p, and is a subsphere of the standard real unit sphere 
S = {u € K" : ||m|| = 1}. This context allows us to study the inverse participation 
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Figure 3. Histograms of IPR values for the Laplacian matrix on random regular 
graphs with n = 1000 vertices and degrees 2 : = 3,4,5,10,20,50. Solid lines 
represent fits to a Gaussian distribution. 


ratio on a space containing by observing that the function IPR{x) is just a 
polynomial of n variables Xi composed from the sum of fourth order monomials 
IPR{x) = n{xf + • • • + xf^) which maps points x G S{p) C K" to K. We may thus 
investigate the origin of the mean value of « 3 of the IPR distributions shown in Fig.[^ 
by recalling that the first moment (the average value) of a function / : K” —)■ K 
over a region is calculated by integrating the product of f{x) and the probability P{x) 
of choosing x from the region of integration i.e. 


p,]pji = IPR{x) = j IPR{x)P{x)da{p) (11) 

J S{p) 

where da{p) is the measure on S{p) and we have introduced the overline notation to 
indicate the continuous average. Similarly the second moment, p-’jpp, corresponding 
to the variance of the IPR distributions in Fig. can be investigated by averaging the 
eighth order polynomial [IPR{x)]‘^ on S{p) and employing the usual identity 


d'lPR = — {dipn) ■ 


( 12 ) 


While there is numerical evidence that the components of an individual 
eigenvector of a random regular graph are Gaussian distributed we proceed with 
the assumption that the additional average over graph realizations that we perform 
will cause P{x) to tend towards the uniform distribution: 


P{x) = 


fs(p) ' 


(13) 


Using Eq. (13) in (11) and (12) thus allows us to compute the continuous first 


and second moments of the inverse participation ratio via the integration of fourth 
and eighth order polynomials on S{p). This can be accomplished using the short note 
of Folland |2y| which provides a formula for integrating a polynomial over a sphere. 
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Folland’s formula is stated for a monomial • • • xfj' and extends linearly to 

polynomials with numerical coefficients Ca- Furthermore, the integral of 

a monomial is dependent only on the Gamma function r(6) = t^~^e~* dt where 

6 is a complex number with positive real part. To proceed, write da for the surface 
measure on the unit sphere S C M”. Folland’s result is the following. 

Theorem 14. Let x‘^ = x^^ ■ ■ ■ xf^ be a monomial, so that aj > 0 for all 
^ ^ j ^ n. Setting bj = ^{aj + 1), 


x^ da = 


if some Oj is odd. 


2r(5i)r(62).-.r(6„) „ 

— 7 "-;-II all Qj are even. 

r(6i + 62 + • ■ ■ + " 


We wish to average the polynomial IPR{x) over the subsphere S{p) C S that 
is orthogonal to the Perron-Frobenius vector p. Thus, in order to use Theorem 
we must first rotate our subsphere S{p) around a (n — 2)-dimensional subspace to 
coincide with the subsphere S{en) as depicted for the case n = 3 in Fig. 



Figure 4. The rotation procedure for n = 3 from an oriented to standard 
subsphere (circle) embedded in 


Note that although the target sphere S'(e„) is defined by coordinates in K", the 
n**' coordinate of each of its points equals zero. Hence, S'(e„) is realized as the unit 
sphere of which is embedded in K" according to the rule (j/i, • ■ ■, 2/n-i) 

(yi, 0). Changing variables therefore allows the direct application of 

Theorem 14 to the unit sphere in to achieve our result. The remainder of 

this section is devoted to these analytic calculations. 


4.I. The rotation matrix Q 

The required change of variables is performed via a rotation matrix Q € SO(n) which 
has the property that for any y S S{en) there exists x S S{p) such that 


y = Qx. 


( 15 ) 
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We begin the construction of Q by using the Gram-Schmidt process to find two 
orthonormal vectors in the plane defined by e„ and p: 


vi = p 




(16) 


V2 


(p ■ ^n) 

I l^n (p ■ ^n) 


\/n(n - 1) 


(17) 


To align p with e„, we need to perform a rotation by an angle 9 defined by: 


cos 9 = p - Bn 


1 




(18) 


around the plane formed by Vi and V 2 and the identity space spanned by the (n— 2)- 
dimensional complement of the orthonormal basis. Hence, we have 

Q = 1 + sin 0 {v 2 Vi — Vi ^ V 2 ) + (cos 0 — 1) (tii G Ui + U 2 G ^ 2 ) (19) 


where 1 is the n x n identity matrix and G represents the tensor product of vectors. 
Considering vector components: vu = ^l\fn and V 2 i = (nSni — 1)/ \/n{n — 1) we 
express the rotation matrix in a form more useful for performing explicit calculations: 


Qij 


1 - y/n 
n — 1 


(1 ^in ^nj “f “f ^in 


Using Eq. (20) it is therefore straightforward to confirm that: 


n 

(i) '^Qij = and 

i=i 

n 

(a) V a; G K" such that ||a;|| = 1 and a; • p = 0, ?/„ = ^ QnjXj = 0 . 

1=1 


( 20 ) 


4-2. Evaluation of the inverse participation ratio moments 

Folland’s straightforward application of Theorem[l4|shows that the (n—2)-dimensional 
measure of each of the spheres 5'(e„) and S{p) is 27iT”“i)/2/r(^^), whose reciprocal is 
P{x), the probability of uniformly choosing a point from such a sphere. To proceed, 
write da{en) for the (n — 2)-dimensional surface measure of the sphere S'(e„) and 
da{p) for the (n — 2)-dimensional surface measure of the sphere 5'(p). Applying the 
multivariable change of basis formula for integrals to Eq. 0 therefore yields 


IPR{x) 


r Y ( —) 

X(P, 

r (—) r 

y^^^^/Pi?(Q-^y)|J(Q-^)|da(e„). 


( 21 ) 


Note that Q is an orthonormal rotation matrix, so that the Jacobian J{Q~^) = 1. 
As mentioned above, j/n = 0 since y • e„ = 0, so that all our monomials have 
(n — 1) variables. Theorem 14 further guarantees that although the expansions of 
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IPR{Q~^y) and [IPR{Q~^y)]'^ have many monomial terms with odd exponents, 
our moment calculations give non-zero values for only those monomial terms having 
all their exponents even. Combining the structure of Q~^ = Q’’’, the multinomial 
theorem, the values of |Appendix A[ and the expansion results of [Appendix B| we now 
give closed forms for the first and second moments of the inverse participation ratio 
over S{p). 

Applying Eq. (21), the first moment is 
_ f 

=2^^ y„. («"!') 


n —1 


= + yly'i 


i—l \k—1 


kl 


' n—1 


n—1 


EI Eet+Eolol. 


4 r{i±3)^ 




kl 


( 22 ) 


where we have used the notation that a prime on a multiply indexed sum enforces the 
constraint that no equal indices are included, i.e. 






and the monomial averages over the subsphere y'^ and y'^y^ have been computed using 
Theorem with the individual results given in Table |A1| The double and triple 
summations over the components of the rotation matrix Q are evaluated in jAppendixj 
l^and substituting Eqs. (B.3| and (B.4) into Eq. (22) yields: 


3n 


h'lPR — 


(n + l)(n - 1) 


29 -h 30v^ -h 5n 24 9 

(1 + V^)2 n 

(i/n-l)(3i/n-|-5)(n-2) 
n{ln + 1)2 


= 3- 


n — 1 


(23) 


which has the observed limiting value of 3 for n ^ 1. 

The calculation of the second central moment proceeds in a similar fashion by 
using Eq. (12) and applying the above general preliminaries to the square of the inverse 
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IPR^ (x) 

r(^) 

27r("-i)/2 

n—1 

k^l 


+ ytyf 

k/ 

n—1 

+ ytyf 


'S(e„) 


[IPR{Q ^y)]^ dcr(e„) 




n—1 


n 

Qki +y]] QkiQkj 

ij 

n " 

28y^(5®jQ^j + + IGQkiQuQkjQej) 

i=l i,j 

n •n 

35y^(5|jQ^j + {QkiQij + ^^QkiQuQkjQ'ij 

ij 

+ ^^QkiQiiQkjQij) 


i=l 


'‘J2'ytyhn 


k^i,m 


n 

210 y] QlQiQL + E' {(^QlQlQ 


2 

mj 


+ 7‘^Q]^iQliQmiQkjQ^jQ'fnj “t” ^^QkiQ^iQkjQijQ 


mj) 


+ 5] ylyiy^nVp 

k,^,m,p 


ipj 


io5EqLQ«OLQp. + E' {^QlQlQlA 

i=l i,j 

+ '^‘^Q^iQiiQmiQpjQejQmj + ^AQkiQuQ iniQpiQkjQijQmjQpj^ 


(24) 


and combining the results of Table X2 and Appendix B we find 


IPR^ {x) = 9 + 


48 

n + 1 


270 210 

—V7 + - 


(25) 


Subtracting the square of the first moment, we arrive at the final expression for the 
second moment of the inverse participation ratio on S{p) 


2 24n(n — 2)(n — 3) 

h-IPR = 


(n + 5)(n + 3)(n + 1)^ 
which indeed tends to zero as n —>■ oo. 


(26) 


5. Comparison With Exact Diagonalization Results 

Having uncovered the origin of the universal number 3 as lim„_,,oo M/p/j foi' fhe 
mean of the continuous IPR, we now undertake a systematic comparison of exact 
diagonalization results for the inverse participation ratio of the Laplacian on finite sized 
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random regular graphs and the predictions from continuous analysis on subspheres 
embedded in K" as a function of z and n. 

5.1. 1st IPR moment 

The finite size scaling behavior of the first moment of the IPR can be quantified by 
explicitly computing the average of the IPR over all non-Perron-Frobenius eigenvectors 
X ^ and then further averaging this quantity over graph realizations. In particular, 
we define the mode-averaged IPR (first IPR moment) for a given graph to be 

p-^{z,n) = ^—'^IPR{x), (27) 

n — 1 

x£S 


where the overline (• • •) now represents an average over the (n — 1) eigenvectors in £. 
Fig. [^depicts the n dependence of Eq. (27) for all graph degrees considered, where we 
have averaged over Nq = 5000 random regular graphs for n < 5000 and Nq = 1000 
graphs for n > 5000. 



Number of Vertices n 


Figure 5. Graph and mode averaged inverse participation ratio vs. the number 
of vertices n for different graph degrees 2 : (symbols). The solid line shows the 
finite size prediction for obtained by averaging over a sphere. 


The solid line describes the function = 3 — 6/(n-|-l) derived in Eq. (23) 

by averaging the IPR polynomial over the subsphere S{p). There is good agreement 
for n > 1000, seemingly independent of graph degree. The error bars are obtained by 
computing the standard deviation of p~^ over all graphs in the generated set, with 
the largest uncertainties occurring for z = 3. We postpone a discussion of the size and 
z-dependence of graph-to-graph fluctuations until the end of this section. 

We may investigate deviations between p~^ and p\pp by defining a normalized 

residual _ 

(P-^) 


Ai(z,n) = 1 - 


iApr 


(28) 


which is plotted in Pig. (left) as a function of n for different values of z. The residual 
decays with increasing n with the correction being well fit by an empirically determined 
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Figure 6. Normalized deviation of the graph and mode averaged inverse 
participation ratio for the Laplacian on ran dom regular graphs from the sphere 
averaged value of given in Eq. |2 m in the text. The left panel shows 

deviations that scale as 1/n with a single ^-independent fitting parameter, as 
confirmed in the right panel. 


function of the form Ci/n where the fitting constant was found to be Ci = 4.15 ± 0.2. 
Within the accuracy of the fitting procedure, we determine that there exists no offset 
in the residual as n —>■ oo and any z dependence of ci is far below the errorbars, as 
seen in the right panel of Fig. 


5.2. 2nd IPR moment 


Next, we consider the prediction of Eq. (26) by studying the second moment of the 
distribution of IPR values on finite sized random regular graphs: 
averaged over Nq = 5000 unique graphs for n < 5000 and Nq = 1000 for n > 5000. 
The results are shown in Fig[^ where now deviations from the sphere-averaged value 
PiPR (included as a solid line) are observed for all values of n and z considered. 


3 

I 


p; 

Ph 


OJ 

o 

a 

I 



Figure 7. The graph averaged second moment of the mode averaged inverse 
participation ratio distribution for random regular graphs of varying degree vs. 
the number of vertices. The solid line shows the sphere averaged value, 


The degree dependence is the most obvious: the exact diagonalization results 
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are systematically larger than fJ-'jpji for small z, with the discrepancy decreasing as z 
increases. We have not included data for z = 3 in Fig. [^as these points lie mostly 
off the scale and the peculiarities of this degree will be carefully investigated in the 
following subsection. Additionally, due to the logarithmic scale, we have only plotted 
errorbars showing the additive uncertainty across graphs. For z > 4, the standard 
deviation is on the order of the symbol sizes. 

We again define a normalized residual for the second moment: 


A2(z,n) = 1 - 


^(p-l)2-(p-l)2) 

h-^IPR 


and the absolute value |A 2 (z,n)| is plotted in Fig. 


(29) 




Figure 8. The normalized positive residual between the second moment of the 
inverse participation ratio distribution for large random regular graphs and the 
value obtained by averaging over a sphere. The left panel shows that the 2 
dependence of the correction can be fit via a second order polynomial in 1/z 
with an offset that persists in the n ^ 00 limit (right panel). 

Here, the dominant deviations from the first sphere averaged value are 

degree dependent, and they can be described by a function of the form c^+c^/z+c^/z^. 
The values of the fitting constants C 2 depend on n, with the solid line in the left 
panel representing their average values for n > 1000. Different htting functions were 
investigated, including those with non-integer negative powers of z, but the resulting 
second order polynomial in 1/z provided the optimal value of the least square fitting 
value. The right panel of Fig. shows the n dependence of A 2 and it appears that 
it remains non-zero even as n —>■ cx). This is consistent with the fitting parameter Cg 
which is finite within errorbars for all values of z considered. 

5.3. Effects of localized eigenvectors 

We now address the issue of the large graph-to-graph variance observed around the 
hrst and second moments of the IPR distribution for small z which led us to remove 
data points corresponding to z = 3 in Figs.[^and[^ This omitted data is displayed in 
Fig. where the graph averaged first and second moments of the IPR distribution are 
shown as a function of n for z = 3. The displayed errorbars correspond to one standard 
deviation and we observe that the effects are most pronounced for the variance of the 
IPR. Data points consistently fall above the sphere averaged value of and the 

mean value between graphs can vary by as much as 1000%. However, for both the first 
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Figure 9. The first (left) and second (right) moments of the inverse participation 
ratio distribution vs. the number of vertices computed via exact diagonalization 
of 5000 random regular graphs with degree z = 3. The errorbars correspond to 
one standard deviation, and they are significantly larger than those observed for 
z > 3. Solid gray lines are the predicted sphere averaged values of and 

pfpR defined in Eqs. ( |23| l and l |26| l respectively. 


and second moments, the data is consistent with values of and computed 
using continuous analysis. The existence of a single outlying data point corresponding 
to n = 5000 that is much closer to in the right panel of Fig. is suggestive 

that the sample set of unique random regular graphs may not be large enough to 
capture the variation in the eigenvector components amongst graphs as measured by 
the inverse participation ratio. 

To better understand the prevalence of this effect for graphs of small degree, 
we have exactly diagonalized the Laplacian for every one of the Me = 4060 unique 
random regular graphs with n = 16 and z = 3 |30] . Analyzing the eigenvectors and 
computing the inverse participation ratio, we find: 


(p-i) = 2.4 ±0.4 


(30) 


= 0.95 ± 1.2 . 


(31) 


The origin of these sizeable graph-to-graph variations is uncovered in the left panel of 
Fig. 10 which shows a histogram of all IPR values (excluding the Perron-Frobenius 
mode) for the complete graph set plotted against their corresponding eigenvalue. The 
frequency of IPR values is shown on a logarithmic color scale from light to dark, and 
we observe two spikes near e = z and e = z ± 1 with the inverse participation ratio 
ranging up to its maximal value of n/2 = 8. In the right panel of Fig. 10 we show the 
maximum value of the IPR across all n = 16 modes and find at least one value of 8 
for nearly 15% of graphs in the set. 

These graph realizations contain special eigenvectors {v} C £ of the Laplacian, 
with IPR{v) = nl2. More generally, a vector with IPR{v) = n/k has exactly k sites 
non-zero, and therefore is of the form 


Vi = 


(- 1 )«‘ 


(32) 


where i{k) G Pk-i,ik} and qi = ±1 such that ~ 0. Clearly 

this is only possible for k = 2m where m G {1, • • • , [n/2j}, with the total number of 
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Figure 10. A histogram with shading on a logarithmic scale of the inverse 
participation ratio vs. eigenvalue for all 4060 random regular graphs with n = 16 
and 2 = 3 (left). The right panel shows the maximum value of the inverse 
participation ratio across all modes for a given graph plotted as a function of a 
fictitious graph index. By relabeling the graph index we also include the maximum 
value sorted by size. 


such vectors given by the multinomial coefficient 


m, n — 2mJ ml ml {n — 2m)! 

For vectors of this form, the inverse participation ratio is indeed given by 


i—l 


(33) 


which is exactly what we observe for k = 2. This maximal value for the IPR 
occurs for the most localized mode that is still compatible with orthonormality to 
the Perron-Frobenius eigenvector and consists of exactly two non-zero values with 
opposite sign. We have confirmed that such eigenvectors indeed appear in our large-n 
graph ensembles for z = 3. Such vectors have a maximal nodal domain count of unity 
[T5] and we observe that almost all eigenvectors with IPR{v) = n/2 have non-zero 
components of opposite sign situated in vector components with consecutive indices. 

The large graph-to-graph variations displayed in Fig. can thus be traced back 
to these localized eigenvectors in combination with the prefactor of n in the definition 
of the IPR given in Eq. By averaging over the sphere, we found in Eq. (23) 


that ~ 3 for n ^ 1. However as just demonstrated, localized eigenvectors can 

contribute IPR values of 0(n) to the first moment. This implies that the variance 
around the mean could contain dominant terms scaling like which will always have 
an effect when averaging over a finite number of large graphs. 


6. Discussion 

In this paper we have investigated the first and second moments of the distribution of 
the inverse participation ratio for all eigenvectors of the discrete Laplacian on finite 
size random regular graphs. By exactly diagonalizing large ensembles of graphs of up 
to n = 10000 vertices we have shown that the average inverse participation ratio is 
well described by the analytic form = 3 — 6/(n -I- 1) obtained by averaging a 
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fourth order polynomial corresponding to the IPR over the sphere S{p) with uniform 
probability. This result is independent of the degree z and we empirically determine 
corrections for small n that fall off as 1/n. Thus, for n ^ 1, we find that the 
first moment of the inverse participation ratio approaches a constant of order unity 
(lim„_>oo(p“^) = 3) for all values of z, implying that the average eigenvector of the 
Laplacian on random regular graphs is delocalized, in agreement with the predictions 
of Ref. [23] for regular graphs with z ~ n. It is interesting to note that a factor of 
three also appears when averaging over the entire domain of the normal distribution, 
which agrees with the observations of Elon m that the components of eigenvectors 
on regular graphs are Gaussian distributed. 

For the variance of the inverse participation ratio over all modes, we have 
again compared our exact random regular graph eigenvectors with an analytical 
result from continuous averaging over S{p) with uniform measure where we find 
h-^iPR — 24n(n —2)(n —3)/[(n + 5)(n + 3)(n+l)^]. Here we observe weaker agreement 
that now strongly depends on the graph degree. This discrepancy appears to persist 
in the limits z,n —)■ oo. When computing the standard deviation of our results over 
an ensemble of up to 5000 random regular graphs, we find that for small values of 
the graph degree z, large fluctuations between graph eigenvectors can cause variations 
in the first and second moments as large as 1000%. By analyzing the complete set 
of graphs for z = 3 and n = 16 we have shown that such deviations arise from 
graphs where the Laplacian has localized eigenvectors consisting of only a few non¬ 
zero elements. Although we have no proof that these vectors appear as eigenvectors of 
the Laplacian for finite size random regular graphs in general, we have demonstrated 
that there are factorially many such eigenvectors that are orthonormal to the Perron- 
Frobenius mode. 

For small values of z, the number of random regular graphs is certainly large, 
but the sample of graphs that we use for computations appears to be large enough 
to capture a significant amount of the variation in graph structure. As mentioned 
above, we have indeed found random regular graphs which have eigenvectors with the 
IPR value of n/2. Since our sample size of graphs is large but fixed, and we conduct 
this fixed size sampling for various values of z and n, the size of the sample pales in 
comparison to the actual number of random regular graphs when n is large. Hence, 
the fact that the errorbars in the mean and variance of the inverse participation ratio 
for large n and z appear to be small is likely due to our samples of random regular 
graphs not being fully representative of the variation which exists in the complete 
set. As z and n increase, extremely large ensembles of graphs need to be studied to 
capture the effects of localized modes, with implications for studying physical models 
with observables computed on regular graphs. 

It would be interesting to explore this issue further, although considerable 
computational resources would have to be employed to diagonalize such large numbers 
of graph Laplacians for n 3> 1. Determining the combinatorial, physical, and 
theoretical significance of these localized eigenvectors is thus left as a topic of future 
work. 
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Appendix A. Averages of Monomials on the Sphere 

This appendix contains results of averaging monomials over S{en) necessary for the 
integral calculations of Section]^ Note that in our case, the polynomial IPR{Q~^y) 
consists of monomials in the variables yi,...,yn-i since yn = 0. Theorem |14| 
guarantees that the integral of a monomial is non-zero precisely when its variables 
have even exponents, and as such, we give values in only this case. 

When performing the first moment calculations we apply Theorem to two 
distinct monomial types. The monomial yf has Oi = 4 for a fixed i and aj = 0 for all 
j ^ i, while the monomial yfy^ has Ui = aj = 2 and ak = 0 ior k i,j. In the first 
case, we therefore have 6^ = 5/2 for a single i and bj = 1/2 ior j ^ i. In the second 
case, bi = bj = 3/2 and 6^ = 1/2 for k ^ i,j. The results of averaging over «S'(e„) are 
gathered in Table [AT] 


Average 

vt 

vlvl 

Value 



r(r^) 

r(r^) 


Table Al. The average values of the degree four monomials with even exponents 
in yi ,..., y-n—i taken over the domain S{en)- 


A similar analysis of the five monomial types appearing in the second moment’s 
polynomial [IPR{Q~^y)]^ yields non-zero averages for those monomials of total degree 


eight which we list in Table A2 


Average 

vl 





vlvl 

ytve 

ylvivli 

ylyhl^yl 

Value 

105 

16 r(^^) 

15 

16 r(^^) 

16 

3(v^r' 

16 r(^^) 

16 rm^) 


Table A2. The average values of the degree eight monomials with even 
exponents in yi ,..., yn—i taken over the domain S{en)- 

We now discuss each of the denominators appearing in these monomial integral 
calculations, and do so taking the value into account. This quantity is the 

numerator of the probability P{x) of choosing points uniformly from the sphere and 
is a prefactor of each moment calculation. A closed form for the first and second IPR 
moments therefore depends on our ability to simplify several n-dependent ratios. 
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For the first moment 
denominator of Theorem 

r(n-i 


2-)/r(^)- 


14 


the total degree of the monomials is four, giving the 
a value of r(^^). Hence, we seek a closed form for 
For a positive integer k, the Gamma function takes the values 
F(fc) = [k — 1)! and F(fc + ^) = (fc — |) • • • so that we have the following 

derivations: 

If n — 1 = 2fc is even, then 


F(^^) {k-iy. 1 4 

F(^^) “ {k + l)\ ~ ~ (n+l)(n- 1)' 


(A.l) 


On the other hand, if n — 1 is odd, then n — 1 = 2fc + 1 for some k and 


r(^) ^ ^ 4 ^ 4 

F(=^) (A:+f)(fc+4)...(4)v/if (2fc + 3)(2fc + l) (n + l)(n-l)- 


The second moment calculati on c ontains only monomials of total degree eight, so 
that the denominator of Theorem [l^ is In a fashion similar to the derivation 

above, we give an alternate form lor the fraction F(^^^)/F(^^^). 

When n — 1 = 2A: is even, we have 


r(^) (fc-i)! _1_ 

r(i^) (fc + 3)! (fc + 3)(A: + 2)(fc + l)(fc) 

On the other hand, for odd n — 1 = 2A: + 1, 


16 

(n + 5)(n + 3)(n + l)(n — 1)' 

(A.3) 


r(^) (fc-i)(fc-i)---(5)v^ 

r(^) + i) (^ + I) (^ + I) ■ ■ ■ (1) 


_16_ 

(n + 5)(n + 3)(n + l)(n — 1) 

(A.4) 


Appendix B. Evaluation of Q Summations 


In this appendix we provide details on the evaluation of the summations over the 


components of the rotation matrix Q given in Eq. (20) that appear in the expressions 
for the first (Eq. (22)) and second (Eq. p4| )) moments of the of the inverse participation 
ratio. These evaluations are performed by first noting that all such powers of Q only 
appear with the first index smaller than n and in this restricted case we can write: 


{Q^3r = 


(n + i/n)' 


[(- 1 )® + as(n)Snj + /3s(n)S^j] 


(B.l) 


where as(n) and /3s (n) are power dependent functions of n that are listed in Table B1 
and we have used the fact that Sfj = Sij and j = 0 Vs, s' > 1 since i < n. 
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s 

0 ^( 74 ) 

Ps{n) 

1 

-^/n 

n + ,Jn 

2 

74 -I- 2^/n 

(74 — 1)(74 -1- 2y/n) 

3 

—y/n (3 -1- 3-^74 -1- 74 ) 

(74 -1- — SyBi —2n + 2rr‘^’^ P r?') 

4 

(?1 H- 2y^)(2 + 2y^ + Tl) 

(74 — 1)(74 -1- 2\fn) 

X (r? -I- 2r?l'^ — 74 — 2\fn -|- 2) 

6 

(74 -|- 2y/n)(l + \/n + 74 ) 

X (3 -I- 3\A4 + n) 

(74 -|- 2\fn)(n — 1)(1 — \/n P 2rC'l’^ P n^) 

X (3 — 3-\/r4 - 274 -1- 274^/^ -|- n^) 

8 

{n + 2\/n)(2 + 2^/n + n) 

X (2 -1- 4-v/n -1- 674 -h 474^/^ -1- 74 ^) 

(n + 2\/n){n — 1) 

X (2 — 2\/n — n + 2 t?^'^ P n^) 

X (2 — 4^/n + 2n + — 5n^ 

- 874 ®/^ -h 2n^ + 474 ”^/^ -h 74"*) 


Table Bl. Expressions for as{n) and /3s{n) appearing as coefficients of 
Kronecker 5-functions when evaluating powers of the components of the rotation 
matrix {QijY where i < n in Eq. ||B.1[|. 


Appendix B.l. 1st IPR moment 


We begin by using Eq. (B.l I to perform the double and triple summations appearing 
in the expression for the first moment of the inverse participation ratio in Eq. (22). 

n n—1 


Qki “ /-t4 “4 + P4in)Ski\ 

i-1 k=l t ' V 7 

[n(n — 1 ) + a 4 {n){n — 1 ) + /34{n){n — 1 )] 


n n— 1 


1 


(n -h i/n)4 
n — 1 
(n + 


n + 04 ( 74 ) -I- P 4 {n)] 


and using the values of 04 ( 74 ) and P 4 {n) in Table Bl we find 

n n— 1 




i—1 k—1 


29 + 30y/n + 5n 24 9 

(1 -I- v^)2 ^ 74 ' 


(B.2) 


(B.3) 
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Now, dropping the explicit n dependence of the as and /3s functions for simplicity, the 
triple summation may be performed in a similar manner: 


n 

i—l k,i 


kiQli 


1 


(n 


i=l 


1 

n—1 

V 

(n 

+ V^)4 

k^l 


1 

n—1 

V 

(n 

+ 

k^^ 

(n 

- ^){n- 

-2) 


EE (1 + a2Sni + /32lJfci) (1 + Oi2^ni + 

M 

n— 1 

(1 + a2)^ + + /32<5fci)(l + h^ti) 


(n + i/n)^ 

(i/n- l)(3A/n + 5)(n - 2) 

n{y/n + 1)2 


- [(1 + 0 : 2 )^ + n - 1 + 2/32] 


(B.4) 


where we have used a 2 (n) and ^(n) from Table B1 


Appendix B.2. 2nd IPR moment 

There are seventeen individual summations appearing in the expression for the average 
of the square of the inverse participation ratio over the sphere given in Eq. (24) and 
we will include the details of only a representative sample here. All can be performed 
using similar techniques employing Eq. (B.l) and Table B1 and we begin with the sum 
over the components of Q® which can be evaluated in exact analogy with Eq. (B.2|: 

n—1 n 




1 


k—1i—l 


(n + ^/ri)^ 
\fn - 1 


EE [1 + CXs^ni Ps^ki 


k—1i—l 


(49 + 7v^ - 7n + 119n®/2 + 2171^ - ISSn®/^ + 


(1 + 

+ llln^/2 - 57n®/2 - 13n® + ISn"/^ + 7n® + . 

The next novel term includes Q® which appears as the third summation and is 
evaluated in a similar method as in Eq. (B.4) albeit with the modification that it 
involves the product of two different powers of Q: 

'fT'—l n ^ n—1 n 

(1 + a^Sni + fdedki) (1 + a25ni + /32^ti) 


1 


(n + v^)® 


n—1 

E 

k^i L 

(n — l)(n — 2) 

(n + i/n)® 

(i/n - 1) (-2 + n) 
(1 + V^)®n3 


n — 1 


(1 + a6)(l + 0 : 2 ) + ^~~^(1 + /36^fei)(l + P25li) 

i=l 

[(1 + ae){l + a 2 ) + n - I + Pa + (32] 

(37 + 3lVn + lOn + 40n3/2 + 2977^ - ISn®/^ 
-14n® + 4n^/2 + 5n^ + . 
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The final type of term contains mixed second indices between the different rotation 
matrix powers and we consider the thirteenth sum in Eq. (241 as a representative of 
this set. The strategy is the same for all such terms and involves performing the 
inner summation by extracting the terms with i = n and j = n and performing the 
summations over j and i, then breaking the remaining restricted sum over i ^ j < n— 1 
into the difference of an unrestricted sum over all values of i,j < n — 1 and one with 
i = j < n — We have 


. 

imj 


{n + y/nf ^ QhQnQkjQijQl 

n—1 

= (~1 + O 3 ) (~1 + Oil) ( — 1 + (flSkj) (—1 + PlSej) (1 + 

i=i 

n—1 

+ (~1 + 0 ^ 1 ) (1 + 0 ^ 2 ) (~1 + (~1 + 

n—1 n—1 

+ (~1 + ^sSki) (—1 + PlSii) ( — 1 + PlSkj) (—1 + PiSlj) (1 + 


2=1 

n—1 




~ (~1 + Ps^ki) (~1 + (—1 +/3l^fei) (1 +/32'Jmi) 


i=l 


— [(~1 + ai)(~l + — 1 — 2/3i + ( 32 ) + (1 + q;2)(~1 + 0!i)^{n — 1 — /3i — P 3 ) 

+ (n — 1 — /3i — I3^){n — 1 — 2/3i + (32) — (n — 1 — 3/3i + /32 — f33 3- (3i(32 + (3f)] 


and putting everything together: 

n — 1 n 


(1 - v^)(n - 2)(n - 3) 


E' E' QlQuQk,Q,,Ql, = - - — - (29 + 397^ + 4n - 28n^/^ 


k,£,m i,j 


(1 + i/n) 


- I2n^ + 12n®/2 + lOn^ + 271^/^) . 
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